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Abstract 

We compute the effective action and correlators of the Polyakov loop operator in 
the Schwinger model at finite temperature and discuss the realization of the discrete 
symmetries that occur there. We show that, due to nonlocal effects of massless 
fermions in two spacetime dimensions, the discrete symmetry which governs the 
screening of charges is spontaneously broken even in an effective one-dimensional 
model, when the volume is infinite. In this limit, the thermal state of the Schwinger 
model screens an arbitrary external charge; consequently the model is in the 
deconfined phase, with the charge of the deconfined fermions completely screened. In 
a finite volume we show that the Schwinger model is always confining. 
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1 Introduction 



The Schwinger model [|^ is an exactly solvable model which describes the electrody- 
namics of a massless fermion in 1+1 spacetime dimensions. It is the classic example 
of a confining gauge theory and displays some of the features of quantum chromo- 
dynamics or other higher dimensional gauge theories where strong infrared effects 
are important 0]. In one space dimension the tree level electron-positron potential, 
V{x) (X e'^\x\, is already confining without quantum fluctuations. Detailed analysis 
in both the path integral [|, |^ and operator methods f^, ||, shows that the spec- 
trum is completely gapped, exhibits chiral symmetry breaking and has no charged 
excitations. 

Over the past few decades, there has been considerable interest in the properties 
of the Schwinger model at finite temperature, both in the massless [H, ^, |10|, |ll|, |l^, 
0, 0, |l^, |16|, |l^, |18|, |l^, [2^] and in the massive I^T] case. The most interesting 



question as to whether the chiral symmetry breaking seen in the vacuum is restored 
at high temperatures was answered in the negative long ago The breaking of 

the chiral symmetry is a consequence of the axial anomaly, rather than spontaneous 
symmetry breaking, and the axial anomaly, being a result of short distance physics, 
is insensitive to temperature. Thus, many of the features of the Schwinger model are 
not changed by temperature. 

In this paper, we advocate the use of temperature to explore the spectrum of 
the Schwinger model. The thermal state is a density matrix with non-vanishing 
contributions from all states in the spectrum with finite energy and thus contains 
information about all of the states. 

As an infrared regularization, we shall consider the space as a circle of circumfer- 
ence L on which all of the basic fields of the Schwinger model have periodic boundary 
conditions. The Hamiltonian quantization in this regularization and some of the ques- 
tions concerning topology and theta-vacua which arise in this case were discussed by 
Manton p3[ . These issues, as well as questions of Wilson line phase dynamics and 



correlators, were further developed in Ref. ||24 



Part of our motivation for this work is to test a recent idea 251 that the 



Polyakov loop operator, introduced by Polyakov [|26| and Susskind [S^] as an order 
parameter for confinement in non-abelian Yang-Mills theory in higher dimensions 
is also a useful operator for Abelian gauge theory. In non-Abelian gauge theory, 
the Polyakov loop has the limitation that it is an order parameter for confinement 
only in models where all of the fields are invariant under global gauge tranforms in 
the center of the gauge group, i.e. are in the adjoint or some other zero 'N-ality' 
representation. In electrodynamics, on the other hand, it can be used in any model 
which is essentially compact in the sense that all of the charges of the dynamical 
fields are integer multiples of some basic charge [25]. Then, as was argued in |25], 
the Polyakov loop with an incommensurate charge can be used to probe the ability 
of the electrodynamic system to screen external charges. 
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1.1 Zn symmetry of finite temperature Yang-Mills theory 

We shall first review the role of the Polyakov loop operator as an order parameter for 
confinement in Yang-Mills theory at finite temperature. This is conventionally seen 
in the Euclidean path integral formulation of the finite temperature gauge theory. In 
that formulation, the Polyakov loop operator measures the holonomy of the gauge 
connection in the periodic Euclidean time, 



P{x) = trPexp li J (irAo(r, x) 



whose correlators in a finite temperature Yang-Mills theory are defined by the 
Euclidean path integral 



{Pix,)...Pixjp^{m)---PKyn)) 

_ J dA,e- 1 P{x,) . . . P{xjp^{m) ■ ■ ■ PKyn) 



(2) 



} dA^e Jo ' 
where the gauge field has periodic boundary conditions, 

A^(l/T,x) =A^(0,x) (3) 

Since the Yang- Mills field transforms in the adjoint representation of the gauge group, 

^) = S-^^l^, ^)9{t. + i9'^{T. ^)V^(r, x)g{T, x) (4) 

they remain periodic under gauge transformations which are periodic up to an 
element, Z, of the center of the group, 

g{\lT,x) = Zg{^,x) (5) 

The center of SU{N) is Z^, the additive group of the integers modulo A^, whereas 
the center of U{N) is the Abelian group U{1). The coset of the group of all gauge 
transformations modulo those which are strictly periodic is a global transformation 
by elements in the center of the gauge group. 

As well as pure Yang-Mills theory, any gauge theory which has matter fields which 
transform in the adjoint, or any other zero 'N-ality' representation of the gauge 
group will have this symmetry of the path integral. Furthermore, this symmetry 
exists for any gauge group which has a non-trivial center. In the following, we shall 
take the gauge group to be SU(N), with the center Z^. The non-trivial topological 
structure arises from the non-vanishing first homotopy IIi{SU{N)/Zn) = Z^. An 
important question is whether or not this Zn symmetry is spontaneously broken at 
finite temperature. 

Under the gauge transformation (|^), the Polyakov loop operator transforms as 

P\x) = ZP{x) (6) 
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Therefore, this operator can be used as an order parameter for breaking of the 
symmetry. 

The connection between the breaking of Zn symmetry and confinement is through 
the fact that the correlators of Polyakov loop operators 

^-Fi,u■■■,^^,y^,■■■,y„)/T ^ ^p(^^) P{xjP\y,) . . . P\y^)) (7) 

can be interpreted as giving the free energy F{xi, . . . ,Xm,yi, ■ ■ ■ ,yn) of the finite 
temperature gauge theory with an array of classical, external, fundamental repre- 
sentation quark sources at positions Xi, . . . ,Xm and anti-quark sources at positions 
yi, . . . ,yn- The normalization of the correlator subtracts the free energy of the gauge 
theory at the same temperature in the absence of sources. 

If the Zn symmetry is not spontaneously broken, the correlator in (|^) vanishes 
unless m = n modulo A^, i.e. unless the quarks and anti-quarks occur in the right 
numbers to make up mesons, which are quark-antiquark pairs, or baryons or anti- 
baryons, which are groups of quarks or anti-quarks, respectively. The vanishing 
of the correlator is interpreted as the quark charge distribution having infinite free 
energy when it has quantum numbers which cannot be combined into color singlets, 
i.e. as quark confinement. On the other hand, if the Z^ symmetry is spontaneously 
broken, the correlators can be non-zero and even quark charge distributions which 
cannot form color singlets can have finite free energy. 

Furthermore, the free energy F{xi, . . . ,Xm,yi, ■ ■ ■ ,yn) can be viewed as the 
effective potential energy of the array of quarks and antiquarks. For example, the 
effective interaction between a quark and anti-quark is given by F{x,y). If this 
correlator increases with distance, as it would in a confined phase where there is a 
non-zero string tension, then the correlator of Polyakov loop operators would have 
the clustering property 

lim e-^(^^'«^)/^ = lim {P{x)P\y)) = (8) 

>oo \x—y\~^OD 

This implies that the Zn symmetry is unbroken and is consistent with the vanishing 
of the expectation value of the loop operator, 

(P(f)) = (9) 

On the other hand, in the deconfined phase, one would expect the quark-antiquark 
potential to fall off to zero with some screening length (the non-Abelian analog of 
Debye screening). In that case, 

lim e-^(^^'«^)/^ = lim {P{x)P\y)) ^ (10) 

\x—y\~*oo \x—y\~*oo 

This implies that the Zjv symmetry is spontaneously broken and is consistent with 
the loop operator having a non-vanishing expectation value 

(P(f)) ^ (11) 
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To properly compute this one-point function, one should as usual introduce a 
symmetry breaking external field through a term in the action such as 

Sx = X Jdx (P(f) + P\x)) (12) 

and compute the one-point function, which would be non-zero when A is not zero. 
Then, the occurence of symmetry breaking would be seen when the limit 

lmi(P(x))A (13) 

is non-zero. 

This formalism is well developed for finite temperature Yang-Mills theory and 
some other pure gauge theories such as compact U{1) and some Z^- gauge theory. 
All non-trivial pure gauge theories in spacetime dimensions greater than two exhibit 
a high-temperature de-confined phase and almost all have a phase transition to 
a confined phase at some critical temperature. Details are summarized in the 



comprehensive review by Svetitsky, pS 



However, for gauge theories with dynamical quarks, the Polyakov loop operator 
is not a useful order parameter to characterize a confining phase. The reason is that, 
since the quark fields transform in the fundamental representation of the gauge group, 
their action and boundary conditions in the path integral are invariant under only 
strictly periodic gauge transformations. Thus, fundamental representation quarks 
are said to break the Zj^r symmetry explicitly. The free energy of a distribution of 
external quarks is always finite. This is interpreted as the possibility of pair creation 
of dynamical quark-antiquark pairs so that the dynamical quarks can screen the color 
of any external distribution of quark or anti-quark sources. 

1.2 Z symmetry of quantum electrodynamics 

In an Abelian gauge field theory, the Polyakov loop operator is defined by the analog 
of (0) 



P{x) = exp li J drAoi^T, 



x) (14) 



and its correlators are computed using the finite temperature path integral of 
quantum electrodynamics 

with the (anti-)periodic boundary conditions, 

A^(l/T,f) =A^(0,x) (16) 
z^(l/T, x) = -^(1/T, x) , V^(l/r, x) = -^(0, x) (17) 
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These boundary conditions, as well as the measure and action in the functional 
integral are invariant under the gauge transformation 



(18) 



(19) 



when the gauge function is periodic up to an additive constant 



X(l/T, x) = x(0, x) + 27rn/ e 



(20) 



where n is an integer. The coset of the group of all allowed gauge transformations 
modulo the group of all strictly periodic gauge transformations is ni(f/(l)) = Z, the 
additive group of the integers. Note that this is a symmetry of the functional integral 
representation of the partition function even in the presence of dynamical electrons. 
It is an interesting and well-defined question to ask whether this symmetry is realized 
in a spontaneously broken or unbroken phase in quantum electrodynamics. 

The Abelian Polyakov loop operator transforms under a gauge transformation as 



and is not invariant under the Z symmetry unless e is an integer multiple of the 
electron charge e. This transformation law was noted by Hansson, Nielsen and Zahed 



Polyakov loop operator transforms under a Z^ subgroup of Z. Thus, as in finite 
temperature Yang-Mills theory, the Abelian version of the Polyakov loop operator 
can be used as an order parameter for the Z symmetry. 

The Z symmetry is related to charge screening and confinement in quantum 
electrodynamics in the same way as the Z^v symmetry of finite temperature Yang- 
Mills theory. The correlators of Polyakov loop operators measure the free energy of 
the electro dynamic system in the presence of a distribution of static charged sources. 
The two-point function, for example. 



measures the effective interaction potential between particles with charges e and 
— e and positions x and y, respectively. In a de-confined phase, we would expect 
Debye screening and the asymptotic form of the potential at large separations to 
decay exponentially with the Debye mass of the photon. This would imply that the 
correlator of two Polyakov loop operators approaches a constant at large separations. 
This implies spontaneous breaking of the Z symmetry. 

In a confined phase, there should be a string tension, and Fg-eix, y) increases 
with separation. This would give a decay of the two-point correlator of Polyakov 
loops consistent with a Z symmetric phase. 



PL{x) = e 




(21) 



1T9| | when the incommensurate charge e was a fraction of the electron charge and the 
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In 3+1-dimensional quantum electrodynamics, we would expect that, at least 
in the physically observed weak coupling regime, the Z symmetry is broken 



spontaneously at all temperatures. It has recently been argued p5| that in 2+1- 
dimensional parity invariant electrodynamics, at least if the electron mass is large 
enough, both the confined and de-confined phases should exist with a Kosterlitz- 
Thouless type of phase transition between them at some finite temperature. In the 
following we shall examine the case of 1+1-dimensional electrodynamics. There, when 
the mass of the electron is non-zero, the dimension of the space is too low to allow a 
phase transition. Based on the results of ||2^ one can expect that the theory exists 
in a confined phase at any temperature (except for some specific ^-vacua found in 

H). 

On the other hand, we will be able to show that the Z symmetry is spontaneously 
broken in 1+1-dimensional electrodynamics when the electron is massless, i.e. in the 
Schwinger model. 



1.3 Z symmetry of the Schwinger model 

In this paper we shall examine the expectation value of the Polyakov loop operator 
in the Schwinger model. The Z symmetry transforms the temporal component of the 
gauge field as 

v4o ^ Ao + 27rnT/e , n^Z . (23) 

We shall give an interpretation for the Z symmetry in the Hamiltonian formalism 
in terms of the quantization of charge, in states of the thermal ensemble. If Z is a 
good symmetry, all charged states through which the thermal system fiuctuates have 
charges which are quantized in units of the electron charge. If Z is spontaneously 
broken, there are quantum states available which have arbitrary charge. If Z were 
broken to a subgroup, Zjq, this would imply that there were fractionally charged 
states with charges quantized in units oi e/N where e is the charge of the dynamical 
electron. An explicit realization of the latter breaking pattern may be of relevance 
for applications to one dimensional condensed matter systems. 

Our analysis of the finite temperature Schwinger model with one fiavor of fermions 
leads us to the following results: 

i. )In the one-dimensional Coulomb gas, which can be regarded as a certain limit of 
electrodynamics which has very massive charged particles, the Z symmetry breaking 
problem resembles that of the quantum pendulum problem, or 1-dimensional sine- 
gordon theory. The Z symmetry is unbroken at all temperatures, corresponding to a 
confining state. 

ii. ) In the Schwinger model where the space is a circle with circumference L and with 
periodic boundary conditions for both the photon and electron fields, we compute 
the expectation value of the Polyakov loop operator and its correlators. We find that 
the expectation value of the Polyakov Loop operator with electric charge e an integer 
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multiple of the electron charge is a non-zero computable constant. When the charge 
in the loop operator is not an integer multiple of the electron charge, the expectation 
value of the Polyakov loop operator vanishes, 

^exp |ie J (iry4o(r, x)| ^ = if e 7^ integer -e (24) 

at all temperatures T. This can be seen as the consequence of the discrete symmetry 
( PBD which is realized in an unbroken phase when the volume is finite. 

iii.) In the infinite volume limit, the Z symmetry is spontaneously broken. This is 
seen by examining the following limits: 

lim f lim / e''^/o'^^'^"^°("'") e-''^/o A constant ^ (25) 

\x—y\~*oo \L— +00 \ / / 

for any charge e and at all temperatures. We interpret this as implying that the 
thermal state of the Schwinger model can screen arbitrary external charges. The 
exact form of the correlator is known and the asymptotic, exponential decay of the 
correlation function is governed by the Schwinger mass of the photon, m| = e^/vr. 

1.4 Symmetry breaking in one dimension? 

The result that the Z symmetry is spontaneously broken was anticipated by Hansson, 
Nielsen and Zahed [113. It is surprising in the sense that, as we shall argue in 



the following sections, the effective action for the Polyakov loop operator is an one 
dimensional field theory with a discrete Z symmetry. Normally such symmetries 
cannot be spontaneously broken, as the long-range correlations described by (|25| ) 
are forbidden by the accompanying strong infrared effects. From another viewpoint, 
the ordered state of the broken symmetry theory is unstable to the condensation of 
domain walls. 

This can be understood by a simple argument: If we consider a one dimensional 

system with N sites and n domain walls, the entropy of the state can be estimated 

f N\ 

by noting that the domain walls can be arranged in I I ways, leading to entropy 



^ = ln (26) 




If the domain wall has energy e the free energy at temperature T for large N and n 
is then given by 

F = ne + f (nlnn + (A^ - n) \n{N -n)-N\nN^ (27) 

Note that for all values of the domain wall energy e and temperature T, the entropy 
always grows faster than the energy as n is increased. This leads to a condensation of 
domain walls. The equilibrium number of domain walls has a Fermi-Dirac distribution 
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If the size of the system is Na where a is the lattice spacing, the correlation length 
is of order the mean distance between domain walls, 

^^{1+ e^/^) a (29) 

which is always small, of order the "lattice spacing" or inverse ultraviolet cutoff. Thus, 
domain wall condensation would seem to always destroy one-dimensional order. 

We shall argue that the Schwinger model evades domain wall condensation at all 
finite temperatures by having domain walls with infinite energy. This occurs because 
the domain walls are actually instantons in a static gauge. The fermion determinant 
vanishes on instanton configurations, giving the instantons an infinite free energy. 
Thus, the only way out of the above argument, that e/T — oo, is actually realized in 
the Schwinger model. 

When the electrons have a mass, one can expect that the domain wall energy 
for small mass diverges logarithmically, e ^ — Tln(m/yu), for small m where is a 
dimensional constant related to the fermion mass m and the confining scale which is 
given by the electric charge e. Thus, if the electron in the Schwinger model had non- 
zero mass the domain walls would have finite energy, the correlation length would 
be 

^^{l + fi/m)a (30) 

and the domain wall condensation would ruin the symmetry breaking at all temper- 
atures, apparently even in the zero temperature limit. 



1.5 Deconfinement versus superconductivity 

There is another interpretation of the physical state of the Schwinger model 
alternative to confinement. The fact that the photon has a mass can be interpreted 
as the Schwinger model being a superconductor or, since in one dimension there is no 
possibility of magnetic fields and therefore no Meissner effect, a perfect conductor. 
This is seen by considering the current induced in the Schwinger model ground state 
by an external electric field which can be obtained from the exact identities for current 
conservation 

'^,{j,{x))a = (31) 
and the axial anomaly equation which can be presented as 

e^.V^{j.)A = l^e^^V^^:^* (32) 

which makes use of the kinematical identity relating the axial and vector currents of 
two dimensional fermions 

(fa) A = ie^,v{ju)A (33) 
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The above equations have the solution 



0,(0:))^ = '-'-^E^-^ (34) 



in terms of the external electric field, This is a superconducting response. 

For example, if the electric field is spatially constant, it has the solution 



(Jo) A = , {Ji)a = -E'^'^'t (35) 

TT 

where the current increases linearly with time. 

This superconducting response can lead to a super-screening of electric fields which 
would otherwise be caused by external charges. We propose this as an alternative 
to the other obvious interpretation of the breaking of the Z symmetry, the loss of 
confinement. 

Our results regarding the Z symmetry breaking support the conclusions of Refs. 
pi] , pD| , p2[ , where the concept of screening versus confinement in 1+1 dimensional 



field theories was discussed. There, deconfinement was interpreted as arising from 
liberation of so called bleached states, with the charges of the deconfined fermions 
being completely screened. The "bleached states" were originally found in the 
operatorial solution of the Schwinger model in [^]. In later works P, ^ they were 
argued to be unphysical in that they are created by non-local operators and therefore 
are unattainable by the action of operators within the algebra of local observables. 
The strong infrared effects driving the theory in infinite volume are essentially non- 
local, however, so we shall interpret our results as indicating the possible emergence 
of "bleached states" in the finite temperature infinite volume limit. 

In the next Section, we shall present two simple examples where the realization of 
the Z symmetry is in the unbroken phase. In the subsequent Section we shall review 
the Hamiltonian formulation of the Schwinger model. It is somewhat independent of 
the rest of the paper and is intended mainly to fix notation and remind the reader 
of the standard picture. In Section 4, we describe the path integral representation 
of the partition function of the Schwinger model at finite temperature. We also 
introduce the effective action for the Polyakov loop operator and make explicit the 
physical interpretation of the Z symmetry. In Section 5, we calculate the Polyakov 
loop expectation values and prove results ii) and iii). in Section 6, we present a 
discussion of our conclusions. 



2 Two Simple Examples 

Before we solve for the Polyakov loop correlator in the Schwinger model, let us 
consider the following examples. 



10 



2.1 Free electrodynamics in 1+1 dimensions 

First, let us consider the case of two dimensional pure U(l) gauge theory. The 
correlator of Polyakov loop operators is given by 



^-Q^.e,// drA^ir,.,)^ ^ JdA,{x)e [[^ 6 ^^^^ 



where the finite temperature path integral is done with periodic boundary conditions. 
The partition function has the formal symmetry 

Aq{x,t) ^ Ao{x,t) + constant (37) 

which, because of the absence of the dynamical electron field, is larger than the Z 
symmetry. 

It is straightforward to perform the gaussian integral in (|36D to obtain the 
exponential of the 1-dimensional Coulomb energy. The result has an infrared 
simgularity unless 

E = (38) 
If this constraint is obeyed, we obtain the expression 



(39) 



This is the usual confining 1-dimensional coulomb potential. It corresponds to a state 
where the symmetry under translation of Aq in the path integral is unbroken. This is 
seen, for example, in the correlator (^e*'^^"^°*-'^'^-*e~*'^^'^''*-'^'^'') which has the asymptotic 
form 

lim ^e^e/Ao{r,x)g-ie/Ao(r,|y)^ ^ jj^ ^_e2|x-y|^Q ^^q^ 

\x~y\—*oo \x~y\-^oo 

The cluster decomposition implies that the symmetry is unbroken at any finite 
temperature. 

2.2 1-Dimensional Coulomb gas 

The 1-dimensional Coulomb gas is implicitely solvable through the Gibbs ensemble 
calculation of Ref. |^^. Here, we obtain the grand canonical partition function of a 



neutral Coulomb gas by the following construction. Consider the statistical mechanics 
of a state with m classical particles with charge e occupying positions xi, . . .Xm and 
n classical particles with charge — e occupying positions yi,...,yn. The partition 
function is given by 

^ (41) 
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We multiply by the statistical factor for identical particles, l/m\n\ and a fugacity 
paramter A™""*"", average over positions Xi and yi by integrating and sum over m and 
n to obtain the partition function 

rl/T 

Z[X, T] = (42) 

If we fix the gauge 

^Ao(r,x) = (43) 
we can do the integral over A^ and obtain the one-dimensional sine-gordon theory 

Z\\T\ = I JldAo{x)exp j dx {^{dAQ{x) / dxf - A cos(eAo(x)/T)^^ (44) 

The effective action for Aq{x) explicitly has the symmetry under the shift 
Aq Aq + 27™T/e In the one-dimensional system (^4]) this symmetry cannot be 
spontaneously broken for any values (aside from zero or infinity) of the parameters A 
and T. Thus, the expectation value of the Polyakov loop must vanish unless it has 
charge e. This we interpret as confinement. There is no confinement-deconfinement 
transition in this model. 



3 Hamiltonian Formalism 

3.1 Hamiltonian, gauge constraints and theta-states 

We shall consider 1+1-dimensional electrodynamics defined on a compact space, 
X G [0,L]. We begin by working in the canonical, Hamiltonian formalism. The 
Hamiltonian is derived from the action 

S = J dt£ dx - ■ (iV + eA)tlj) (45) 

In this action, the canonical momentum conjugate to the spatial component of the 
gauge field {Ai{t,x), which we shall shortly rename A{t,x)) is the electric field 
E{t, x) = -Foi(t, x) = VtAi — VxAo- The momentum conjugate to the fermion ip(t, x) 
is iip'^it, x). The non- vanishing equal time (anti-) commutation relations are therefore 

[A{x),E{y)] = i5{x~y) 
[^{x),i,\y)} = 5{x-y) (46) 

The temporal component of the gauge field, Aoit^x), appears in the action (|45|) 
without time derivatives and acts as a Lagrange multiplier to impose the constraint 
of invariance under gauge transformations. The Hamiltonian is obtained as 

H = £ dx (^^E\x) + ip\x)a{'iV + eA{x))i;{x)^ (47) 
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where a = 75 = 70 71 is a 2 x 2 Hermitean Dirac matrix and V = d/dx. The massless 
Dirac Hamihonian can be decomposed into Hamiltonians for left and right movers as 

i^Dirac = dx (^i(x)(zV + eA{x))^L{x) - ^lj(a;)(^V + eA{x))^R{x)) (48) 
All fields have periodic boundary conditions in space, 

A{L) = A{0) E{L) = E{0) 

tpiL) = V(0) -^HL) = iP^O) (49) 

The Hamiltonian and commutation relations must be supplemented by the first class 
constraint, or Gauss' law, which is the operator obtained by taking a functional 
derivative of the action (|45| ) by Ao, 

G{x) = -VE{x) - ei)\x)i){x) ~ (50) 

and ensuring that the quantum states are invariant under time-independent gauge 
transformations. The latter are generated by the operator 

G[x]= rdxx{x)G{x) (51) 
Jo 

where x{x) is a periodic function, x{L) = x(0)- The action of the operator (|5T|) is 

A{x) + Vxix) 
E{x) 

e'^''^''^ij{x) 

e--xW^t(a;), (52) 

This is a symmetry of the Hamiltonian ( ^7D and of the commutation relations (^Bj), 
which preserves the boundary conditions (^). 

There is a larger class of gauge transformations under which the Hamiltonian 
and commutation relations are invariant, which preserve the boundary conditions 
(^91) and which are not generated by the Gauss operator G[x]- These have gauge 
functions which are not strictly periodic, but obey the condition 

Xn{x + L) = Xn{x) + 27rn/e . (53) 

This guarantees that both the electron operator and the gauge field boundary 
condition is unchanged. Such 'large' gauge transformations can always be expressed 
as a periodic gauge transformation plus a representative of the large, non-periodic 
transformations as 

Xnix) = Xoix) + 2Txnx/Le , (54) 

where Xo{x) is periodic. 



e^^MA(x) e-*^M = 

e^^W E{x) e-^^M = 

e^^MV'(x) e-^^M = 

e^^M V^Vx) e-*^M = 
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Large gauge transformations are implemented by an unitary operator exp{iGe[xn]), 
where the Hermitean operator is 

Gi[Xn] = £ {Vxn{x)E{x) - exn{x)ij\x)^{x)) (55) 

Using (|5^), this generator can be written as a large gauge transformation generator 
and a Gauss' operator 

GAXn] = G[xo] + -r- / dxE{x) - —- / dxx^\x)ij{x) (56) 
Le Jo L Jo 



The quantization of the model with commutator algebra (|i6|), Hamiltonian (|4y| ) 



and constraint (BUp can proceed in two different ways. First, one can solve the 



constraint ( pOf ) at the classical level by imposing an auxiliary condition on the 
remaining degrees of freedom. The second and equivalent approach, which we shall 
pursue in the following, is to quantize the dynamical system specified by (|^) and 
(^6]) as it is. Then, on the Hilbert space which represents the algebra (|46|) and where 
the Hamiltonian is diagonalizable, we shall impose the physical state condition 

G[xo] I physical state > = (57) 

The physical states are thus invariant under all periodic gauge transformations. 

However, they need not be invariant under the set of all gauge transformations. 
In fact, it is only necessary that they transform under a unitary irreducible 
representation of the coset of time-independent gauge transformations modulo the 
periodic ones. The coset group is isomorphic to the translation group of the integers, 
Z, whose unitary irreducible representations are one dimensional phases, e*^" . Thus, 
if we implement a large gauge transformation using the operator the physical 

states should transform as 

QiGdxr.] Iphysical state, e> = e'""^ |physical state, 6 > (58) 

In this way, the physical states are characterized by a theta-angle. 

Like the theta-angle of non-Abelian gauge theories in four spacetime dimensions 
P3| there exists a canonical transformation which removes the theta-angle from 



the states and introduces a theta term in the action. The unitary operator which 
implements this transformation is 



U{e) = exp (^-iee/27r £ dxA{ 



x) (59) 



In the new system, the theta angle is absent from the physical states and the electric 
field operator is altered. The new hamiltonian is 

H = £ dx(^\{E{x) + ee/2-Kf + i)\x)a{iV + eA)^{x)) (60) 
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In this way, one sees that theta has the interpretation of a constant background 
electric field, as discussed by Coleman, Jackiw and Susskind I^Qj. In the original 



action, taking into account the role of A^, (|45|) is modified by the addition of the 
conventional 6 term, 

S = J (fx (iF^, + ^Foi - ■ («V + eA)ij{x)) (61) 

It turns out that, in the massless Schwinger model, the physical states do not 
depend on 6. It is also possible to see that the parameter 6, which appears with 
the topological term, in the action is invariant. In the following we shall retain the 
theta-dependence in order to demonstrate the theta-independence of the partition 
function. 

In the next Section we shall discuss the construction of the path integral 
representation of the thermodynamic partition function. 

4 Path integral representation of the partition 
function 

The thermodynamic description of field theory is most conveniently obtained from the 
partition function which for a gauge theory is gotten by taking a trace over physical 
states of the Gibbs distribution, e~^^^ , where T is the temperature and we work 
in units where the Boltzmann's constant as well as the Planck's constant and the 
speed of light are equal to one. In constructing the partition function it is convenient 
to consider all the states which represent the commutator algebra (|46|) and insert a 
projection operator which projects over the physical states, and onto a sector with a 
fixed vacuum angle 6. The trace is thus given using a complete set of states which 
represent (^), 

Z[T] = Y.<s\e-''/^Pe\s> (62) 

s 

The appropriate projection operator is constructed from the unitary operator which 
generates gauge transformations 

^-.—Ee-Zfeje-W (63) 

where we have integrated over all time independent gauge transformations and 
divided by the (infinite) volume of the gauge group. This results in the expression 
for the partition function 

Z[T] = ^ E e-^"^ / [dXn] e-^-[>^"] (64) 

where the effective action for the gauge degrees of freedom is given by the twisted 
trace 

g-SeffM = ^ < s| e-^/^ > (65) 
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This effective action fias a standard patfi integral representation; in pfiase space 
it is given by 

e-SefiM = / n Yl dil^{T,x)d^\T,x)dA{T,x)dE{T,x) e"^^!'^''^^'^'-^! (66) 
j;e[o,L] Te[o,i/r] 

wfiere tfie Euclidean first order action is 

Se = dr dx [iEA + ^E^ - [V^ + iaV^ + eaA] V'j . (67) 

The electric field E{t,x) has open boundary conditions in time, and the other 
integration variables have twisted (anti-) periodic boundary conditions, 

A{l/T,x) = A{0,x)-Vxn{x) 
ip{l/T,x) = - e^"^"(^V(0,a;) 

V'^(l/r,x) = V^(0,x) e-^^^"(^) (68) 

The Gaussian integral over the canonical momentum E{x) is performed to yield, up 
to a temperature and volume dependent but Xn-independent overall factor, 

e-5eff[xn] ^ f Yl II diP{T,x)diP\r,x)dA{r,x) e"^^!^'^^'^] (69) 

x£[0,L]Te[0,l/T] 

where now the Euclidean action is given by 

Se = j^'^ dr £ dx (i(i)2 - ip^ [V^ + iaV^ + eaA] ^) (70) 

The boundary conditions can be untwisted by a non-periodic gauge transforma- 
tion. This is what normally restores Ao, the temporal component of the gauge field to 
the Euclidean path integral. A suitable non-periodic gauge transformation redefines 
the integration variables as 

A{r,x) ^ A{T,x)-Vxn{x)TT 

^{t,x) ^ e-^"^"(^)^^V'(T,x) (71) 
jjj\r,x) ^ V'^(t,x) e'^>^"(^)^^ 

Note that the spatial boundary conditions for the fermi fields are now changed. The 
resulting path integral now has the action 

rl/T rL 

Se = dr dx [\{Af + \T^{yXnf - [V^ - ieTxn + ^aV^ + eaA] 



(72) 



and the boundary conditions are 



i^il/T,x) = -m^) Hr,L) = e2™^-^(r,0) 

iIj\1/T,x) = ~ilj^{0,x) ij\r,L) = c'^^"'^^ ij\r,0) (73) 

A{l/T,x) = A{0,x) A{t,L) = A{t,0) 

Xn{L) = Xn(0) + 27rn/e 
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In order to compute the effective action for Xn{x) we must compute the path 
integral ( |66D with the Euchdean action (^) and the boundary conditions (ffSl). 
The effective electromagnetic field tensor is given by 

Foi (r, x) = A{t, x) - VAo{t, x) (74) 

where we identify the temporal component of the gauge field in a static gauge 
{VrAo = 0) as 

Ao{x) = Txn{x) (75) 
This field configuration has instanton number n, as it can be seen from the integral 



e 



2tx 



l/T pL 



JO 



— dr dxFQi{T,x) =n (76) 



where we have made use of the fact that the field A{t,x) has periodic boundary 
conditions in both r and x. Thus, the effective vector potential fields in the ?7,'th 
sector are just those which are n instantons in a static gauge. 

4.1 Z symmetry 



We consider the following change of integration variable in the path integral (^, : 



ij{T,x) ^ e2"'=^"V'(r,a;) 

^\t,x) ^ e-2"''=^"V'^(r,x) . (77) 

When k is an integer, the boundary conditions ([73|) are unchanged by this substitution 
and the Jacobian in the path integral measure is one. The net effect is to replace the 
variable Xn{x) by Xn{x) + 2nk/e. Thus, the effective action for Xn has the symmetry^] 

ScslXn] = Scs[xn + 2Txk/e] (78) 

This is a large gauge symmetry analogous to (|53|) which is associated with the periodic 
nature of the space in the temporal rather than spatial direction. However, being 
associated with Euclidean time, it cannot be a basic symmetry of the theory, it 
is rather an effective symmetry of the Euclidean path integral. We shall presently 
discuss its interpretation in the Hamiltonian formalism. 

In order to obtain a physical interpretation of this symmetry, we consider a 
modification of electrodynamics where there is an array of static external charges 
ei,...,efc located at positions Xi,...,Xk- This can be taken into account by a 
modification of the Gauss' law to 

k 

- VE{x) - ei)\x)i){x) - ej5{x - Xj) ~ (79) 



^Note that this could in principle be only a formal symmetry of the path integral. Here, it survives 
path integration because of anomaly cancellation, similar to the cancellation of gauge anomalies. 
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The sole effect of this modification in the partition function ( pi)) is the replacement 

k 



e """^^"j — >■ e 



n e-^"^-^"(^^) (80) 



This implies that the k-point correlator 




(81) 



is the ratio of the partition function for the electrodynamic system in the presence 
of the external charges to the partition function of the same system in the absence 
of external charges. Thus, the free energy of the system with charges, compared to 
that without is given by 

F{ej,x,)-Fo = -Tln^n e'^^^^^^^'^^ ^ (82) 

where the bracket ( ) is the average over the fields Xn with the measure 
En 6xp(— S'cfr[xn] + in9) . Thus, the correlators of the exponential operators mea- 
sure the Coulomb energy of external charges. In this way they probe the ability of 
the system to screen external charges. 

The symmetry of the effective field theory under the Z transformation, if it is not 
spontaneously broken, poses a restriction on the correlators which can be non-zero - 
and therefore it restricts which arrays of external charges can have finite free energy 
(|82D. In finite volume, this symmetry is certainly realized canonically and the result 
is that any expectation value of the form (|8TD averages to zero when the charges do 
not add to multiples of the electron charge, 

^ ei = me . (83) 

i 

Whether this symmetry persists in the infinite volume limit L ^ oo, is an interesting 
question which we shall discuss in following sections. 

From the definition ( ^5]) of the generators of large gauge transformations, we see 
that the Z transformation (|78D changes the generators by 

Gl [Xn] ^ Gi [xn] - 27ck dx . (84) 

JO 

Accordingly, going back to the definition of the partition function in (|6^), we see that 
the Z transformation corresponds to the replacement of the density matrix e~^^^ 
by the operator 

e"^/'^ exp \-27rik dx ^^^] (85) 



in the trace over the physical states. Since all of the physical states obey Gauss' 
law, with finite L and periodic boundary conditions, they have zero fermion number. 
Therefore, the exponential of the fermion number is trivially the unit matrix on the 
space of physical states. The question arises whether this fact persists in the infinite 
volume limit, or if there are states with arbitrary fermion number. 
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4.2 Computing the Fermion Determinant 

Now we want to calculate the effective action for x, the time-component of the gauge 
field in a static gauge ([75|) . For this we first integrate out the fermions from the path 



integral ( |66| , [72D . 

When the fermion mass is zero, the Dirac operator has zero modes for any of the 
field configurations with n 7^ 0. Thus, all the terms in the sum over n except the 
term with n = vanish. 

The Dirac operator has zero modes as a consequence of the Atiyah-Singer index 
theorem (see e.g. Ref. ||3^). An explicit demonstration in the Scwhinger model can 
be found in Sachs and Wipf [jT3|. The Dirac operator has exactly \n\ zero modes in 
the n-instanton sector. 

For expectation values of Polyakov lines, we need to take into account only the 
zero instanton sector, as there are no fermions in our correlators to soak up the zero 
modes. The fermion determinant in the zero instanton sector isQ 

' d^jdiP e/ ^T(v-*eA)v> = det (7 ■ V - ie7 ■ A) (86) 



We begin with a Hodge decomposition of the gauge field 

= + e^^V^f] + 2'Kh^/e (87) 

where are the harmonic modes, 

e T /"VT rL 

K = 7rT / d^M-r^^) ■ (88) 

/TT L Jo Jo 

As the gauge field lies in the zero instanton sector, the pure gauge (0) and coexact 
(Q) part have to be strictly periodic in space, c.f. (ffSl). 

Using this decomposition, the fermion determinant can be written as 

det (7 ■ (V - ieA)) = det ( e^^+^^^ 7 ■ (V - 27Tih) 6"^'^+^^^ ) 

= det (7 ■ (V - 2mh)) det e^^'^^ (89) 

Here, we have assumed that the determinant of e**^ is the inverse of the determinant 
of e""'^ . This can be shown to be true using a gauge invariant regular izat ion, e.g. 
Pauli-Villars regularization. The other factor is 

det ( e^^^^ ) = exp iti2-f5n) (90) 

This is a standard computation; the coexact part of the gauge field carries the chiral 
anomaly, which can be integrated using any gauge invariant regularization. Noting 
that Afl = Fqi, we get 

det(e2^5^) = exp(^^ J d^xn{x)Foi{x)^ 

= exp f / rf^Foi^Foi ) (91) 



27r 

■^Note that compared to the Minkowskian 7-niatrices used in (^), the 7:3 used henceforth obey 
Chfford relations with an EucUdean metric, thus absorbing the extra i in Action ([72|). 
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In the field strength Foi there is no contribution from the harmonic modes of A^. 

The part of the determinant of the Dirac operator containing only the harmonic 
components of the gauge field has the form 

det (7 ■ (V - i27rh)) = J] [((2n + l)vrr + 27ihof + (27rm/L + 27ihif] (92) 

m.n 

which is the square modulus of the chiral determinant 

det D+ = Y[ {{2n + l)nT + 2TTho + 2TTim/L + 27ri/ii) (93) 

mn 

The latter has the Z-symmetry invariance ho ^ ho + kT, k & Z, and depends only 
on the complex coordinate u = ho + ihi. One can compute the determinant ( PHj ) by 
means of a regularization that preserves the Z-symmetry but breaks the holomorphic 
factorization of ( P^ ) PH| . Namely, one shall obtain for ( P^D a Z-invariant result that 
will depend also on the coordinate u. Such a result is the well-known expression of 
l3f) in terms of Jacobi theta and Dedekind eta functions, 



1 



e 



hn 



(0,^t) (94) 



r]{it) ^ J - hi 

where the modular parameter is it = iLT, and the harmonic modes are rescaled with 
the corresponding length scales to get dimensionless quantities: 

h^ = ho/T ; hi = Lhi . (95) 

For the theta functions, we follow the labelling conventions of Mumford , 



e 



a 
b 



Q{z, it) = J2 (-^^^^ + 2mnz^ (96) 
{z,it) = e-"*"'+2"'"("+^)e(2 + 6 + ita, it) . (97) 



As announced the /ig = (u + 'u)^/4 term in (l9^ ) breaks the holomorphic factorization 
of the determinant (p2|). This is a finite local counterterm that can be added to the 
effective action in order to obtain a gauge invariant result [RS]. As a matter of fact, 
the gauge symmetry on the harmonic component ho of the gauge field is nothing but 
the Z-symmetry. Alternatively, one can maintain the holomorphic factorization of 
the determinant (^2|) and loose the Z-symmetry invariance . Our choice is to keep 
the Z-symmetry invariance. 

The total fermion determinant is obtained by combining the coexact piece 
with the harmonic piece, which is the modulus square of the chiral determinant (| 
Expressing the field strength in terms of the spatial gauge field A and the static 
temporal field we finally get for the effective action 



e 
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dhi 



hi 



(0,2t) 



m 



(9J 



X / [dA] exp <j I / dxdr M f 1 - ^ 



The hat on the fields means that the harmonic part is removed. We have denoted 
the Scwhinger mass by m| = e^/vr. 

The gauge choice ( |75| ) has not fixed completely the gauge, since the harmonic 
part of the gauge field is unaffected by (|7^). Consequently the Z-symmetry is still 
present as a residual gauge invariance. 

The hi and A integrations can now be done, yielding (up to normalization) 



^exp|-r/2 dx{{Vxf + ml{x + 27^/6)'} 



(99) 



In order to get a finite result, one has to gauge fix the residual (spatial) large 
gauge transformations by restricting the integration over hi to the period [0, 1] of 
the integrand. 

Equation (|99|) provides the form of the effective action which explicitly realizes 
the symmetry (ffSl). 



5 Polyakov loop correlators 

Now we are in position to calculate expectation values and correlators of Polyakov 
loops, and accordingly to decide, whether the Z symmetry (^) is spontaneously 
broken. As indicated in Section |4.1| , inserting Polyakov loops probes the response of 
the theory to static external charges. We take the external charges to have charge e. 

Due to the zero modes of the Dirac operator, only the zero instanton sector 
contributes to the Polyakov loop correlators (|8TD . The expectation value divides in a 
global and local part: 

^giex(x)^ = 1/Ary[rfx] e-^°ffM e^^'^(^) (100) 
The global part is easily calculated by Poisson resummation:'^ 

dh ^ Q~2nt{n+ho)^+2Trihoe/e ^ \ — [ dh ^ ^-■K/2tu'^+2TviK{e/e+v) 



* Note that as a gauge fixing of large spatial gauge transformations, the domain of the integration 
has to be restricted to the period of the integrand, which depends on the value of e/e. 
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2t ^ ,11 e/e t ^-j^Q^^ 

[ , otherwise 

The local part can be expressed as 



where we used the harmonic oscillator Green's function on the circle (with the 
contribution of the zero- modes subtracted), 

1 Q2TTin(x~y)/L cosh. i\Lm s{l — 2 \x — y\ / L) 



L ^^m1 + At^'^tt? / 2ms sinh Lm| 

(103) 

Combining ( p.01| ) and ( |102| ), we get for the Polyakov loop expectation value 

^ eie-x(x) ) = I exp {-4^ coth iLms} , if e/e e Z ^^^^^ 
[ , otherwise 

This is consistent with the expected unbroken realization of the Z symmetry for 
finite volume. The heat bath screens only external charges that are integer multiples 
of the electron charge, by bounding dynamical fermions from the heat bath with the 
external charge. This can be viewed as a proof of the confining nature of finite volume 
one dimensional electrodynamics. 

In infinite volume, one would expect the contribution of the harmonic modes to 
be trivial. Accordingly, one would expect the nonvanishing value of the Polyakov 
loop in Equation (|104|) to be valid irrespective of the value of e, which would be a 
signal of Z symmetry breaking. 

In order properly to investigate the possible symmetry breaking in the infinite 
volume limit, we need the Polyakov- anti- Polyakov loop correlator. This is again 
readily calculated using the effective action (PP|). The harmonic contributions cancel 
between the Polyakov and anti-Polyakov loops, and we are left with the integration 



= exp[-B'^/T {K{x,x) - K{x,y))} (105) 
f fcoth^Im cosh [iLmsjl- 2 \x-y\/L) y 

In the infinite volume limit we get 
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This result shows that the system does not cluster decompose: 

^ ^iexix) ^-iexiv) ^ l^^^- |-^^| (107) 

This is consistent with the value of the infinite volume Polyakov loop expectation 
value anticipated from ( [104| ). For e = e, our results ( |104| , |106|) agree with the results 



of on Wilson loop expectation values of the zero temperature Scwinger model on 
a circle, upon a Wick-rotation. 

Equation ( |106|) shows that, in the infinite volume limit, there is off diagonal long- 
range order, and the Z symmetry (^) is thus spontaneously broken. The thermal 
state of the Schwinger model can screen arbitrary external charges, and it is in a 
deconfined phase. 

We interpret our results as indicating the possible presence of the "bleached" 
states of Ref. in the spectrum of the Schwinger model. These are states where the 
charges of the deconfined fermions are completely screened by the thermal state. Our 



results thus support the screening vs. confinement discussion of Refs. ||3T], |30|, p2| . 

In Equation (pS]) we interpreted the Z symmetry as imposing a condition on the 
allowed fermion numbers of the states in the theory. As the symmetry is broken in 
infinite volume, we conclude that in infinite volume states with arbitrary fermion 
number may exist. The screening state, being able to screen any charge, not only 
multiples of the fundamental charge, is a non-local, polarized vacuum like state, and 
as such it can a priori have any fermion number. 

6 Concluding Remarks 

In this paper we compute explicitly the effective action and the correlators of the 
Polyakov loop operator in the one flavor Schwinger model at flnite temperature in 
order to investigate the phases of one dimensional Q.E.D. Our aim is to provide a 
convincing test of the recent proposal p5| that the Polyakov loop operator is indeed 



useful to distinguish between a confined and a deconfined phase of an abelian model 
coupled with fermionic matter. 

We present a form of the finite temperature effective action which explicitly 
realizes the Z symmetry. We show that in one-dimensional Q.E.D. with massless 
fermions the Z symmetry is not broken in finite volume. The Z-symmetry is broken 
— due to strong infrared effects — only in the infinite volume limit where there 
is off-diagonal long range order and the physical states may have arbitrary fermion 
number. 

In infinite volume, the thermal state of the Schwinger model can screen an 
arbitrary external charge and, therefore, it is in the deconfined phase. Our explicit 
computation of the Z symmetry breaking in the Schwinger model is supported by 
two simple arguments providing a sound physical intuition for the breaking of a 



23 



discrete symmetry in a one dimensional field theory. The massive Schwinger model, 
on the other hand, is confining and the Z symmetry is not broken, at least when the 
temperature is much greater than the electron mass and the confinement scale is set 
by the dimensional electron charge. 
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